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The eternal inflation scenario predicts that our observable universe resides inside a single bubble 
embedded in a vast inflating multiverse. We present the first observational tests of eternal inflation, 
performing a search for cosmological signatures of collisions with other bubble universes in cosmic 
microwave background data from the WMAP satellite. We conclude that the WMAP 7-year data 
do not warrant augmenting ACDM with bubble collisions, constraining the average number of 
detectable bubble collisions on the full sky TVs < 1-6 at 68% CL. Data from the Planck satellite can 
be used to more deflnitively test the bubble collision hypothesis. 



Introduction: The inflationary paradigm has been 
very successful at explaining the initial conditions giving 
rise to our observable universe. Considering the initial 
conditions for inflation itself leads to the possibility that 
our observable universe might only be a tiny piece of a 
vast multiverse. In this scenario, known as eternal in- 
flation (for a review, see e.g. Ref. !1 ), our observable 
universe resides inside a single bubble nucleated out of 
a false vacuum de Sitter space. The rate of bubble for- 
mation is outpaced by the accelerated expansion of the 
inflating false vacuum, and therefore inflation does not 
end everywhere. 

Eternal inflation is ubiquitous in theories with extra 
dimensions (string theory being the primary example) 
and positive vacuum energy. However, testing this sce- 
nario is extremely difflcult since eternal inflation is a pre- 
inflationary epoch: any signals from outside of our bub- 
ble would naively appear to be stretched to unobservable 
super-horizon scales. While this is in general true, one 
prospect for probing this epoch lies in the observation 
of the collisions between vacuum bubbles. These colli- 
sions produce inhomogeneities in the inner-bubble cos- 
mology, raising the possibility that their effects are im- 
printed in the cosmic microwave background (CMB) [2]. 
In this paper we describe a robust algorithm designed 
to test the hypothesis that there are bubble collisions in 
the Wilkinson Microwave Anisotropy Probe (WMAP) 7- 
year data [3]. More generally, our analysis pipeline can 
be adapted to test any theory predicting localized signa- 
tures on the CMB sky. 

Recent theoretical work (see the review Ref. ^ and ref- 
erences therein) has established that bubble collisions can 
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produce detectable signals and still be compatible with 
our observed cosmology; in addition, there are models 
that predict an expected number of detectable collisions 
Ns larger than one. The value of Ns is highly dependent 
on: the scalar field potential(s) that drive eternal infla- 
tion (which controls the rate of bubble formation); the 
duration of inflation inside our bubble (the more inflation 
the weaker the signal); and the particular realization of 
the CMB sky and bubble collisions we might observe (i.e. 
even a strong signal could be obscured by foregrounds). 
Although there is ample motivation to consider the eter- 
nal inflation scenario, a concrete model (of, say, the string 
theory landscape) providing all of these details does not 
currently exist. 

Nevertheless, a single bubble collision produces a 
rather generic set of signatures that we target in our anal- 
ysis: 

Azimuthal symmetry: A collision between two bub- 
bles leaves an imprint on the CMB sky that has azimuthal 
symmetry. This is a consequence of the S0(2,l) symme- 
try of the spacetime describing the collision of two vac- 
uum bubbles. 

Causal boundary: The surface of last scattering can 
only be affected inside the future light cone of a collision 
event, which forms a ring on the CMB sky. The observed 
temperature of the CMB need not be continuous across 
this boundary. 

Long-wavelength modulation inside the causal 
boundary: A bubble collision is a pre-inflationary relic. 
The effects of a collision are stretched by inflation, and 
induce an overall modulation of the CMB temperature 
anisotropies [5]. 

At the time of last scattering, the signal is long- 
wavelength, has approximate planar symmetry (from the 
small observed value of the normalized curvature density, 
fife), and has a causal boundary. The observed tempera- 
ture fluctuation (to lowest order in cos 9, centered on the 
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north pole) takes the form 
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where zq is the amphtude of the temperature modulation 
at the centre of the collision, Zcrit is the discontinuity of 
the temperature at the causal boundary 9 = ^crit, and 
0(^crit — 9) is a step function. An example is shown 
in the top left of Fig. [T] A temperature modulation of 
this form, with z^rit = 0, was first derived by Chang 
et. al. [S], who analyzed the behaviour of the inflaton 
field in a background thin-wall collision geometry. Our 
model allows, but does not require, the existence of a 
temperature discontinuity at the edge. Including the lo- 
cation of the collision centre, {9o, ipo}, a collision can be 
described by five parameters: {zq, Zcrit,9crit,(^Oj<t'o}- A 
complete model of eternal inflation would predict a prob- 
ability distribution for the number of collisions and their 
parameters. In the absence of such specific predictions 
we use Ns as a continuous parameter that characterizes 
particular models of eternal inflation. The concordance 
ACDM cosmological model is given by the special case 
in which Ng — 0. 

Data and simulations: We calibrate our analysis 
pipeline and determine its sensitivity using simulations 
of the WMAP 7-year data. WMAP has measured the 
intensity and polarization of the microwave sky in five 
frequency bands. We perform our final analysis on the 
foreground-subtracted 94 GHz W-band WMAP temper- 
ature map, as this combines high resolution (determined 
by the detector's 0.22° beam) with small foreground con- 
tamination. To minimize the effects of the residual fore- 
grounds we cut the sky with the conservative KQ75 mask, 
leaving 70.6% of the sky unmasked. 

The WMAP data-reduction and foreground-removal 
pipelines may leave behind systematic errors that we are 
unable to include in our likelihood function. To cali- 
brate our detection thresholds and quantify the expected 
false detection rate due to such systematic effects, we 
use an end-to-end simulation of the W-band channel of 
the WMAP instrument, provided by the WMAP Science 
Team. This is generated from a simulated time ordered 
data stream (containing Galactic foregrounds, CMB fiuc- 
tuations, realistic instrumental noise, finite beam size, 
and other instrumental effects) which is processed using 
the same pipeline applied to the actual data. 

To determine the sensitivity of our analysis pipeline 
to bubble collisions, we generate simulated temperature 
maps containing CMB fiuctuations, a bubble collision, in- 
strumental noise, and smoothing to emulate the W-band 
beam. We consider collisions with 0crit = 5°, 10°, 25° 
and choose centres in a high- and low-noise unmasked 
region of the sky. For each 6'crit and location, we gener- 
ate 35 simulated collisions with parameters logarithmi- 
cally distributed in the ranges 10"^ < zq < 10~^ and 



— 10~^ < Zcrit < — 10~^; the lower limit of the Zcrit 
parameter range produces templates with edges which 
are undetectable in WMAP data. The response of our 
pipeline depends only on the absolute value of zq and 
Zcrit, so the choice of sign for zq and Zcrit is arbitrary. 
We repeat this process for three realizations of the back- 
ground CMB fiuctuations, yielding a total of 210 simu- 
lated sky maps for each of the three collision sizes. 

Method: Our primary goal is to determine, given the 
WMAP 7-year data, what constraints can be placed on 
Ns and whether models predicting > should be be 
preferred over models predicting Ng — 0. From Bayes' 
theorem (e.g. Ref. [6]), this information is contained in 
the posterior probability distribution for Ng given full-sky 
CMB data. Evaluating this would require the computa- 
tionally intractable tasks of inverting the full-sky full- 
resolution CMB covariance matrix and integrating the 
bubble-collision likelihood over a many-dimensional pa- 
rameter space. However, as we describe in more detail in 
a companion paper [Tj, it is possible to approximate the 
full-sky Bayesian evidence integral by a patch-wise anal- 
ysis, if it is possible to first identify the regions on the sky 
that contribute the most to the integrand. Our analysis 
pipeline automatically locates these regions, avoiding a 
posteriori selection effects [8]. Once a set of candidate 
patches has been identified, it is also possible to apply 
further tests of the data in parallel. 

Our analysis pipeline is schematically depicted in 
Fig. [T] A "blob detection" step using needlets (a type 
of wavelet) is used to identify significant features in the 
temperature map and determine their approximate lo- 
cation and angular size. Two parallel verification steps 
are then performed: first, an edge detection algorithm is 
applied to these features to search for circular temper- 
ature discontinuities (zcrit in the template). Second, a 
pixel-based Bayesian model selection and parameter es- 
timation analysis is performed on the regions of the sky 
highlighted by the needlets. These results are used to 
construct an approximation to the full-sky posterior dis- 
tribution of Ns, from which we can determine if the data 
supports augmenting ACDM with bubble collisions. 

Blob detection: To identify the most promising can- 
didate signals in a temperature map, we perform a suite 
of spherical needlet transforms 
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where j = frequency, k — HEALPix pixel [TIT, are 
the spherical harmonics, a£m are the spherical harmonic 
coefficients of an input map, and &^(j) is a frequency- 
dependent filter function determining the needlet shape. 
Two classes of needlet shapes ( "standard" |9] and "Mex- 
ican" |!Tl]) are used to ensure that we are sensitive to 
a wide variety of modulations. We optimize the needlet 
response to a variety of collision templates by adjusting 
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FIG. 1. The signatures of a bubble collision at various stages 
in our analysis pipeline. A collision (top left) induces a 
temperature modulation in the CMB temperature map (top 
right). The "blob" associated with the collision is identified 
by a large needlet response (bottom left), and the presence 
of an edge is highlighted by a large response from the edge 
detection algorithm (bottom right). In parallel with the edge- 
detection step, we perform a Bayesian parameter estimation 
and model selection analysis. 

the band- width of bi{j). Optimal unbiased maximum- 
likelihood estimators of the a£„iS 12 are used at low £ to 
partially offset the effects of the mask on low-frequency 
needlet coefficients; cut-sky a^mS are used at larger 

In the lower left quadrant of Fig. [T] we plot for fixed 
j at each pixel k obtained for a simulated map contain- 
ing a collision. The needlet coefficients clearly take their 
largest value in the vicinity of the collision. Finding the 
frequency of maximum needlet response yields informa- 
tion about the angular scale of the collision. 

For a purely Gaussian uncut CMB sky, the average 
needlet coefficient is zero for all j, k, and the variance 
of the needlet coefficients for fixed j at each pixel k is 
identical and directly related to the scalar temperature 
power spectrum. Cutting the sky introduces a j- and 
fc-dependent bias, and we determine the significance of a 
needlet coefficient by 

gauss, cut I /Q\ 
^Jk = , ^ = , l-i) 

y (/^jfc) gauss, cut 

where the average (/3jfe) gauss, cut and variance 
(/^jfc)gauss,cut are calculated at each pixel from the 
needlet coefficients of 3000 collision-free Gaussian CMB 
realizations with the WMAP 7-year KQ75 sky cut 
applied. 

Because there are many independent needlet coeffi- 
cients in any given realization, there can be features of 
fairly large significance even in a purely Gaussian map. 
In addition, residual foregrounds and artifacts of experi- 
mental systematics can be picked up by the needlet trans- 
form. To minimize the number of such false detections, 
we use the needlet response to the end-to-end simulation 
to define a set of frequency- and shape-dependent detec- 
tion thresholds (in the range 3.00 < Sjk < 3.75). Regions 



with five or more pixels whose needlet coefficients exceed 
these thresholds, and which do not intersect the Galactic 
mask, are grouped into "blobs", and are passed to the 
following steps. A total of 10 false detections are found 
in the end-to-end simulation using these thresholds. The 
angular scale of these features is read from a lookup ta- 
ble, built by finding the range of ^crit over which each 
needlet frequency yields the maximum response to a set 
of simulated collision templates 

We have tested the performance of the blob detection 
step using the full set of simulated collisions. The results 
are presented in Fig. [2] (left panel) as "exclusion" and 
"sensitivity" regions of the {zq, Zcrit} parameter space for 
&cTit — 10° (the 0crit = 5° and 25° results are nearly 
identical) . If we detect all six simulated bubble collisions 
at a given point {zq, Zcth} of parameter space, then a non- 
detection in the real data would imply we can exclude 
such collisions; if we detect a collision for a subset of 
the simulations, then we are sensitive to collisions in this 
region of parameter space. 

Edge detection: This step of the pipeline tests fea- 
tures highlighted by the blob detection stage for circu- 
lar temperature discontinuities (controlled by Zcrit)- To 
generate a set of candidate edges, we have developed 
an implementation of the Canny algorithm [ 13] for the 
HEALPix pixelization scheme. In this algorithm, the 
gradient of an image is generated, smoothed to reduce 
the effects of pixel noise, and thinned into local maxima, 
the strongest of which are stitched together into edges. 

An adaptation of the Circular Hough Transform 
(CHT) algorithm [Mj is then used to assess whether these 
candidate edges fall on circles of varying centre and angu- 
lar radius. If a map of the "CHT score" (the fraction of 
pixels in an annulus surrounding a given center that are 
candidate edges) is sharply peaked at some angular scale 
6, then there is evidence for a circular edge of this radius 
centered on the peak. In the example in the bottom right 
quadrant of Fig. [T] the peak in the CHT score at the true 
centre of the collision shows an unambiguously detected 
edge. We have verified that none of the false detections 
flagged in the end-to-end simulations have a significant 
peak. 

The sensitivity and exclusion regions from the edge- 
detection step are shown in Fig. [2] (right panel) for 
the 10° simulated collisions. For the 5° collisions, the 
proliferation of degree-sized features in the background 
CMB affects performance, yielding a similar sensitivity 
region, but a much smaller exclusion region (including 
only zq > 10^* and Zcrit ^ — 10~*). The exclusion and 
sensitivity regions for the 25° collisions are, again, nearly 
identical to those for the 10° collisions. 

Based on these results, the first two steps of our 
pipeline can detect bubble collisions with central modula- 
tions of j^ol ^ 3x10^^ or causal edges of |zcrit| ^ 3x10^^ 
with 6'crit > 5°. 

Bayesian analysis: Assuming that the bubble colli- 
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FIG. 2. Exclusion (black) and sensitivity (grey) regions for the needlet (left) and edge-detection (right) steps of the analysis 
pipeline applied to a set of Scrit = 10° simulated bubble collisions. Collisions in the exclusion region would be definitively 
detected as long as they were not significantly masked. Collisions in the sensitivity region could be found if they were in a 
favorable location of the sky (i.e. low noise, or a region with a specific realization of CMB fluctuations which did not obscure 
the causal edge). 



sion model likelihood is peaked in the iVb blobs identified 
by the needlets, the (unnormalized) posterior probability 
for the number of detectable collisions, iVg, given a CMB 
data set d covering a sky fraction /sky, can be approxi- 
mated as 



Pr(iVs|d, Uy) cx Pr(7Vs) 6"/-^^= ^ 



(/skyiVs 
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bi,b2,...,b]sf =1 
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where Pr(A^s) is the prior probability of Ng (assumed to 
be uniform) , the pre- factors reflect the fact that the num- 
ber of collisions present on the observable sky, N^, is the 
realization of a Poisson-like process (of mean /sky-^s), 
and pb = Pr(dh|l)/Pr(dfc|0) is the ratio of the model- 
averaged likelihoods (or evidences) of the bubble model 
and ACDM evaluated within a candidate collision region 
(with data sub-set d^). The posterior can therefore be 
built from local measures of how well the data are de- 
scribed by the standard model with and without a single 
collision template. 

We use the nested sampler MultiNest I ') to calculate 
the single-blob evidence for the bubble collision model, 



Pr(db|l) = y"dmPr(m)Pr(d,,|m,l), 



(5) 



by marginalizing the bubble collision likeli- 
hood Pr(db|m, 1) over the model's n parame- 
ters m. The (exact) pixel-space likelihood is 
Pr(db|m,l) cx exp[-(db - t(m))C,-Hdf, - t(m))T/2], 
where t(m) is the temperature modulation caused 



by the collision, and C;, is the pixel-pixel covariance 
matrix of the collision region (computationally limited 
to patches of radius < 11°), including cosmic variance 
given by the best-fit Ci as well as the W-band noise 
and beam from the WMAP 7-year data release. The 
evidence for ACDM, Pr(dfc|0), can be calculated by a 
single evaluation of the bubble likelihood with t(m) = 0. 

The parameter prior Pr(m) in Eq. [5] is derived from 
theory, previous experiments and, as only detectable bub- 
ble collisions are considered, the limitations of the data- 
set and pipeline. Lacking a detailed theoretical predic- 
tion for the amplitude parameters {zq, Zcrit}, we assume a 
uniform prior over the ranges — 10~* < {zq, Zan} < 10"'*. 
Bubble collisions are expected to be distributed isotropi- 
cally on the CMB sky, so we assume uniform priors on the 
full ranges of {cosOq, (po}. The prior range for detectable 
bubble collisions is restricted by CMB power at small 
scales and computational requirements at large scales. 
Based on this, we assume uniform priors on ^crit values 
between 2° < 6'crit < 11°. 

The evidence ratios found in the end-to-end simula- 
tion are bounded by Inpb < —6.6. Evaluating the pos- 
terior, it is maximized at iVg = 0, and we conclude that 
the residual foregrounds and systematics from the end- 
to-end simulation do not provide false support for the 
collision hypothesis. Evidence ratios for simulated colli- 
sions in the needlet and CHT exclusion regions are very 
large {Inpt ~ 100), yielding a full posterior (even for one 
blob) that clearly favors the bubble collision hypothesis 
for a variety of Ng. For simulated collisions in the sensi- 
tivity region, pt can cover a wide range (—7 ^ Inpb < 9), 
and does not always yield full-sky posteriors that favor 
the bubble collision hypothesis. Thus, while we might 
rule out these features as being due to systematics or 
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foregrounds, better data would be needed to definitively 
establish the bubble collision hypothesis. 

Analysis of the WMAP 7- year data: Applying the 
needlet step of our pipeline to the WMAP 7- year W-band 
temperature map with the KQ75 mask, we find 15 sig- 
nificant features, four of which intersect the main Galac- 
tic mask and are hence discarded (discarding patches 
from the evidence integrand leads to a conservative lower 
bound on the integral). The number of features and 
their significance range is fairly consistent with the end- 
to-end simulation (we find needlet significances in the 
range 3.37 < Sjk < 4.58). The edge-detection results 
are also consistent with those from the end-to-end sim- 
ulation, with no clear peaks in the CHT score for any 
of the significant features. We can therefore rule out 
the presence of any large-angular-scale bubble collisions 
with values of {zoj^crit} in the CHT exclusion region of 
Fig. [2j The evidence ratios for the 11 significant features 
are bounded by Inp < —3.8, yielding a posterior for Ng 
that is maximal at Ag = 0. Therefore, the WMAP 7-year 
data-set does not favor the bubble collision hypothesis for 
any value of Ag; from the posterior we find Ng < 1.6 at 
68% probability. 

Interestingly, however, the Bayesian evidence ratios 
obtained for four of the features in the WMAP data 
are systematically larger than expected from false de- 
tections, as calibrated using a collision-free end-to-end 
simulation of the WMAP experiment. Data from the 
Planck satellite [16], which has a resolution three times 
better than that of WMAP, with an order of magnitude 
greater sensitivity, will greatly improve the pipeline's di- 
agnostic power. As we have shown, the non-detection of 
a bubble collision can be used to place constraints on the- 
ories giving rise to eternal inflation; however, if a bubble 
collision is verified by future data, then we will gain an 
insight not only into our own universe but a multiverse 
beyond. 
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